RESONANCES IN TWISTED QUANTUM WAVEGUIDES 
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Abstract. In this paper we consider embedded eigenvalues of a Schrodinger 
Hamiltonian in a waveguide induced by a symmetric perturbation. It is shown 
that these eigenvalues become unstable and turn into resonances after twisting 
of the waveguide. The perturbative expansion of the resonance width is 
calculated for weakly twisted waveguides and the influence of the twist on 
resonances in a concrete model is discussed in detail. 



1. Introduction 

Quantum and electromagnetic waveguides have been studied since many decades, 
see [TBI [22l [23j [27] and also [24] where the classical and the quantum pictures 
are compared. In this framework the spectral analysis of differential operators 
in tubular domains has become a research field of a certain interest (6] [HI [16]. 
Moreover, with the introduction of nano-devices, such as nanotubes, then new 
open problems in quantum transmission for such structures appeared [7]. 

We consider here a waveguide type domain £1 = K x u> (see Figure [1] on the 
left), where the cross section lu of the waveguide is an open bounded and connected 
subset of R 2 . We impose Dirichlet boundary conditions at the boundary of O. 
The spectrum of the free operator — A in L 2 (Vl) is absolutely continuous and covers 
the half- line \E\,oa), where E\ is the lowest eigenvalue of the Dirichlet Laplacian 
on u). It is a well known fact that this spectrum is unstable against perturbations; 
indeed, a negative perturbation, vanishing at infinity, of —A will induce at least one 
bound state below the threshold E\ and new embedded bound states in the half- line 
[Ei,+oo) (see, e.g., Figure[2]). The perturbation can be either of a potential type 
or of a geometric type, see [HI US] and references there. Similar effects occur also 
in two-dimensional strips, [UIU[T5]- These new bound states below the threshold 
Ei correspond to the particles (electrons) which do not propagate along fi, but 
remain localized in a bounded region of f2. 

Recently it has been shown [14] that the presence of bound states in fl can be, 
up to certain extent, suppressed by another geometrical perturbation: the so called 
twisting which is defined as follows. For a given x € K and s := (y, z) G uj we 
define the mapping 

/ ( :Mxw^l 3 
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by 

f e (x, s) = (x, y cos(e a{x)) + zsin(e a(x)), z cos(e a(x)) — j/sin(e a(x))), (1) 

where e > is a real parameter and a : M — > M is a diffcrentiable function. 
Furthermore, we introduce 

fi e := /«(«)■ 

Clearly, f2 e is a tube which is twisted unless the function a is constant (in Figure [T] 
we plot, respectively, a rectangular tube without and with twisting). The result of 
[H] shows that if the cross section uj is not rotationally symmetric and the tube f2 
is twisted, even only locally, then the bound states for the perturbed Hamiltonian 
—A + V do not appear for any negative potential V(x), but only if V is strong 
enough. In other words, one could say that a twisting of a tubular domain f2 
improves the transport of charged particles in f2 in the sense that it protects the 
particles to get trapped by weak perturbations. The repulsive effect of twisting 
has been recently observed also in [17j , where the absence of discrete eigenvalues in 
tubes which are simultaneously mildly curved and mildly twisted. Moreover, in [3J 
the repulsive effect of twisting is demonstrated for bounded tubes whose thickness 
goes to zero. 

Similar results were also obtained for two-dimensional waveguides with com- 
bined boundary conditions or with a local magnetic field [131 HO]. However, the 
geometrical perturbations of the waveguide generically induce also the existence of 
resonances, i.e. metastable states with very long lifetimes, see [101 1111 [T2"] . 

It is the aim of the present paper to describe the influence of twisting on the 
resonances in the waveguides. More precisely, we study the situation in which the 
free Laplacian is perturbed by an attractive potential V(x), which decays at infinity 
along the waveguide direction, where x € R represents the coordinate along the 
waveguide direction. The point spectrum of the perturbed Hamiltonian — A + V(x) 
consists, in addition to the bound states below E\, of infinitely many eigenvalues 
embedded in the continuum [Ex, oo) (see Figure [2]). It was shown in [To] . 

for two-dimensional waveguides, that these embedded eigenvalues generically 
turn into resonances in the presence of a constant magnetic field. Following the 
method of [10] we show that this happens also when the magnetic field is replaced 
by the twisting, provided the cross section lo is not rotationally symmetric, see 
Theorem 1. For weak twisting we also give the perturbative expansion of the 
corresponding resonance width. 

In order to obtain a precise estimate on the imaginary part of the resonances and, 
in particular, to prove that it is strictly negative we consider in Section[5]a concrete 
model in which the potential V approximates a one-dimensional point interaction. 
For such a model we explicitly calculate the leading term of the imaginary part of 
a chosen resonance, see Proposition Q] and we prove that for suitable values of the 
parameters (see Remark|9|), the imaginary part of the resonance is strictly negative. 

2. Preliminaries 

Throughout the paper we will denote by (-,-)h the scalar product in a Hilbert 
space H with the convention (au,v)n = cx(u,v)h for all a € C and u,v e H. 
For a real-valued measurable bounded function V(x) on R we formally define the 
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Figure 1. On the left, the plot of the surface of a rectangular 
waveguide without twisting; in the right, plot of the surface of the 
twisted rectangular waveguide. Bold line represents the boundary 

of U). 



Hamiltonians 

H° = -A and = - A + V(x) in L 2 (n e ) 

with Dirichlet boundary conditions at 90 e . The operator is associated with 
the closed quadratic form 

Qebl>] ■■= [ [IV^I 2 + V(x)\^\ 2 }dxds, (2) 
Jn e 

with the form domain D(QY) = Wj(fi e ). 

Given a test function tjj £ Co°(fi) it is useful to introduce the following shorthand, 

d T tp := yd z ip - zd y ip. (3) 

As usual in such situations, in order to analyze the operator we pass from 
the twisted tube f2 e to the untwisted tube f2 by means of a simple substitution of 
variables. This gives 

QXW\= I (\V s TP\ 2 + \d x 4> + ea(x)d T ijj\ 2 +V{x)\^\ 2 )dxd Sl 
Jn 

with the form domain D(QY) = Wj(f2) and with the notation 

V s ^ := {d y ip,d z ip) . 

In other words, the operator , associated with QY and unitarily equivalent to 
HY 1 acts on its domain in L 2 (il) in the weak sense as 

H V e = -d 2 - 8 2 - [8 X + ea(x) 8 T ] 2 + V(x) = H% + C/ e y , 

where 

H% = -d 2 x d 2 d 2 + V(x) 
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and 

UY = -[dx + ea(x) d T } 2 + d 2 x 

= —2ea(x)d x d T — ea(x)d T — e 2 a 2 (x)d 2 . 

Remark 1. The term UY is a symmetric operator on L 2 (il) with Dirichlet bound- 
ary conditions at dQ. 

In order to show that the embedded eigenvalues of Hq turn into the resonances once 
the waveguide is twisted, we employ the method of the exterior complex scaling in 
combination with the regular perturbation theory [§]. We start by locating the 
spectrum of the untwisted model. 

3. Spectrum of HY 
We will suppose that V satisfies the following 
Assumption A. The function V(x) is not identically equal to zero and 

[ (I + x 2 ) \V(x)\dx < oo and [ V(x)dx<0. (4) 

JTg. JR 

It then follows from [35] (see, e.g. Theorem XIII. 110 in and its Notes) that the 
operator 

h := -flg + V{x) in L 2 (M) 

possesses finitely many negative eigenvalues {^j}jLi, N > 1, each of multiplicity 
one. We denote by tpj(x) the corresponding normalized eigenfunctions. The essen- 
tial spectrum of h covers the positive half- line [0, oo). On the other hand, it is well 
known that the operator — A 1 ^,, i.e. the Dirichlet Laplacian on to, is positive definite 
and has purely discrete spectrum. Let {En}^^ be the non-decreasing sequence of 
its eigenvalues and let Xn(s) denote the associated normalized eigenfunctions. The 
set of such eigenfunctions is an orthonormal basis of L 2 {u>). We denote by 

S = {E = fjLj +E n , j = 1, . . . , N, n > 1} 

the set of eigenvalues of HY with associated normalized eigenvectors 

^n,j(x,s) = ipj(x)Xn{s), 

and 

S + = En [£i,+oo), £_ = En (-oo,J5i) 

where S_ is not empty since fij < for any j. Then, by the standard arguments, 
[25], the spectrum of 

HY — -A + V(x) , ini 2 (Rxw) 

is given by cr(HY) = crd(HY) U <j ess (HY), where 

<J d (HY) = S_ and a ess (HY) = [E 1 ,oo). 

In addition, HY possesses point spectrum embedded into the continuum given by 
£+ (see Figure^). 

We expect that when e becomes non-zero then these embedded eigenvalues generi- 
cally turn into resonances, which are the main object of our study. 

Remark 2. Since the operator HY commutes with complex conjugation then its 
eigenfunctions -0 can be assumed to be real- valued. 
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Figure 2 . The discrete spectrum of Hq consists of finitely many 
simple eigenvalue below E± (denoted by full circle); the essential 
spectrum is given by the half- line [E\, +oo). Furthermore, a non 
empty set of simple eigenvalues (denoted by empty circle) embed- 
ded in the half-line [E\,+oo) occurs. 



4. Complex scaling 

Henceforth, we will employ the method of exterior complex scaling to the oper- 
ator HY ■ In order to do so, we will need some assumptions on the functions V 
and a: 

Assumption B. V extends to analytic function with respect to x in some sector 

M fj := {C S C : | argCI < /?}, with (3 > 0. 
Moreover, V is uniformly bounded in Mp. 

Assumption C. a extends to analytic function with respect to x in 

Mp = Mp U {C € C : |3ffl < /?} , with [3 > 0. 
and a is uniformly bounded in Mp. In addition a(x) > 0, Vx G K. 

Remark 3. Since a is uniformly bounded in Mp then from the Cauchy theorem 
it follows that a is uniformly bounded in A4p> for any < /3' < j3. 

In analogy with [10] we introduce the mapping Sg, which acts as a complex dilation 
in the longitudinal variable x: 

(S e ip){x )S ) = e e / 2 ^(e e x,s), 6eC. 

The transformed operator then takes the form 

HY(9) = SgHYSg 1 = HY(8) + UY(0) , 

where 

HY(9) = SeHYS^ = -e~™ dl - dl - d\ + V(e e x) 

and 

UY{6) ^SgU^Sg 1 = -2ee- e a{e e x)d x d T -e e~ e &{e 9 x) d T -e 2 a 2 (e e x) d 2 T . (5) 

Lemma 1. Let V satisfy assumptions [21 and then Hq (6) is an analytic family 
of type A with respect to 9. Furthermore, the spectrum of (8) has the form (see 
Figure 0i 

a(HY(e))=\J[E n +e-^R+]. (6) 

n 

More precisely, the essential spectrum of Hq (9) consists of the sequence of the half- 
lines E n + e _2lSe IR + , n = 1,2,..., and the discrete spectrum of Hq (9) consists of 
the set of eigenvalues fij + E n with associated eigenvectors 

[ffinM] = [S<fil>nj] {x,s) = e e / 2 <p 3 (e e x)xn(s). (7) 
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Figure 3. The discrete spectrum of Hq(0) consists of a sequence 
of real and simple eigenvalues (denoted by circle); the essential 
spectrum is given by the half-lines E n + e~ 2lise M. + . 

Proof. It follows from Assumption [B] that the family of operators Hq(9) in analytic 
of type A with respect to 9, see [19l Chap. 7]. For what concerns its spectrum it is 
enough to remark that the operator 

h(9) = SghSg 1 = -e~ 29 d 2 x + V(e 9 x) 

in I/ 2 (M) has the spectrum given by 

a{h{9)) = {/xi, w} U e' 2 ^ e R+ . 

By the Ichinose's Lemma, see [25j Sec. XIII. 10], we then obtain the formula ([6]) for 
the spectrum of the sum h(9) — d 2 — d 2 . □ 

Lemma 2. Let V satisfy assumptions^^ and\B\ and let a satisfy assumvtion [Cl 
then the operator UY (8) is a relatively bounded perturbation of H^{0) . Moreover, 
the family of operators fLY (6) is analytic of type A for all 9 such that \9\ < R e , 
where R £ — ► oo as e — ► 0. 

Proof. In order to prove the Lemma we consider a test function ip S (f2) . Using 
the assumptions on V we find out that there exists a positive constant c (here and 
below, c will denote a positive constant whose value changes from line to line) such 
that 

\\hX(8) i;\\ 2 >c{\\d 2 x M\ 2 + lia.^11 2 + RV|| 2 ) - cii^n 2 . 

From this inequality and using the fact that u> is bounded we arrive at 

\\d 2 ^\\ 2 < c(||9 z V|| 2 + ||^|| 2 + ||9^|| 2 + ||^H 2 ) 

< c(||^|| 2 + ||9 2 ^|| 2 + ||^|| 2 ) 

< C (|KW|| 2 + ||^|| 2 ) 

and 

RV|| 2 < c(||9^|| 2 + H^ll 2 ) < c (\\HY(0) V|| 2 + U\\ 2 ) . (8) 
As for the mixed term in Uj(0) we note that 

Ha^ll 2 = (dM, d x d T ^) L2(n) < \\d 2 ^\\ ||a 2 ^|| < \ [||a 2 ^|| 2 + \\d 2 M 2 \ . (9) 

Collecting all these estimates we finally conclude that there exists a positive con- 
stant C, depending on to and V, such that 

\\U? (9) V|| 2 < C (e + e 2 ) 2 (||<(#) V>|| 2 + || V|| 2 ) . (10) 

To prove the second statement of the Lemma we first notice that by assumption 
Ewe have D(H%{6)) = D(H^(0)). By assumption O and [Ml Sec. 7.2] it thus 
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suffices to show that both d x d T and d T are relatively bounded with respect to 
H^{9). However, this follows from © and (10]). □ 

Lemma [2] tells us that the eigenvalues of H^ (9) are analytic functions of 9. By a 
standard argument, [8J, it turns out, that they are in fact independent of 9. The 
non-real eigenvalues of HY (9), for 9 such that 36* > 0, are identified with the 
resonances of HY [8] . 

Remark 4. As a result of the previous proof it follows that UY (9) is a regular 
perturbation of the operator Hq (9). This enables us to apply the analytic pertur- 
bation theory to the eigenvalues of the operator Hq (9). 

Theorem 1. Let E = E n + fij G E+ be a simple eigenvalue of Hq 9) where 9 is a 
fixed complex number such that Q9 > 0. For any ball B centered in E there exists 
s* > 0, depending on E, such that for any e with \e\ < e* , there is an eigenvalue 
E(e) of HY (6) belonging to B and with the imaginary part given by 

SE(e) = -e 2 a + 0(e 3 ) (11) 

where a is a constant independent of e and equal to 



k<k* 



(drXn,Xk)m^\ (vj,^f(E ~ E k ) Vj } L2(m) . (12) 



Here 



Vj = (—2ad x + &)ipj, k* = max{/c : E k - E < 0} 

and Qt stands for the imaginary part of the reduced resolvent [19] of h = —d% + V 
with respect to the eigenvalue \ij : 

(vj,^f(E- E k )vj) L2(m = 
1 

' ' " " " "' L ' ' ' ' >L 2 ( 



lim 



(vj,r(E- (E k +ip))vj} L2(m - (vj,r(E- {E k -ip))v j ) 1 



Proof. Let ip(9) = ip n j(9) be the associated normalized eigenvector ([7]) belonging 
to E, where n and j are fixed. We apply the regular perturbation theory saying 
that for some fixed r > small enough and for any e with modulus small enough, 
in the given ball B r (E) exists only one eigenvalue E(e) of H^ (9) with associated 
eigenvector 



ne) = ^-.l [c-HY(9)y 1 mdc- 



>dB T 

Furthermore, the regular perturbation theory also yields that this eigenvalue is 
given by means of the convergent perturbative series (see, e.g. )25[ XII. 6]) 

($(9),HY(9U £ (9)) r2ln , °° 

where, as usual, 

en = E and ei = , -, — — — = — ; ; — 

are constant with respect to 9; tp is the real-valued vector for 9 = (see Remark 
[2]). These constants eo and e\ are real- valued since UY is a symmetric operator 
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(see Remark[T|). If we prove that SJe2 = — e 2 a + 0(e 3 ) for some a > independent 
of £ then the stated result follows. To this end we recall that (see [TO] ) 

3e 2 = 9fo a (1 + 0(e)) , 

where 

« 2 = -^<f {m,uY(d)[c-H^(e)]- 1 uY(e)m) 1 



t 2 (0) C - E 



= lim f(6,E + ip) = lim f(8 = 0,E + ip) 

p^0+ p->0+ 

and 

mo = -(m,uy(d)[c-HY(e)]- l uY(9)m) 



L 2 (fi) 



+ 

Hence 



<^),c/ £ y (^)) L2(f2) (C-i?)- 1 



a 2 = -(^, 17/ JJ(B + iO)U* V)l 2 (o) := lim -(V, + %p)UXf)^a) (13) 

where R(Q is the reduced resolvent of Hq with respect to the eigenvalue E, see 
[251 X.II.6] and [13 . Recalling that V has the form 

1p(x, s) = 1p n ,j(x, s) = <Pj(x)Xn(s) 

for some n and j and that {Xfc( s )}fcli is a basis of L 2 (uj), we obtain 

oo 

UYi>(x,s) = ^2d k (x)xk(s), where d fc (x) = (xk , i/j) L 2 {w) . 

k=l 

Using the fact that U^-ip S L 2 (f2), Lemma [2] and the dominated convergence 
theorem we conclude that 



ai = ~ ( d k^(E - E k + iO)d k ) L2(m , (14) 
fe=i 

where r(£) is the reduced resolvent of h = —d 2 + V with respect to pj. Concerning 
the imaginary part of a 2 we point out that only finitely many terms on the r.h.s. of 
(fT4| have a non zero imaginary part. The latter follows from the fact that (d k , r(E— 
Ek + iO)dk) is real for any k large enough; more precisely for any k > k*, where 

k* = max{fc : E k - E < 0} . 

From 

d k (x) = £Vj(x)(d T Xn,Xk)L 2 (u) [1 + 0(e)] , Vj = {-2ad x + a)(fj 
we can thus conclude that 

a 2 = -e 2 A[l + 0(e)} , 

where 

A = A n,j = X! \( d rXn,Xk)L^(Lo)\ 2 (vj,r(E - E k + iO)vj) L2{K} 

k<k< 
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is independent of e. This implies (p~2|) since 



a = $sA = ^2 \{drXn,Xk)L^(oj)\ 2 (vj,$sr(E- E k )vj) 



L 2 (R) ' 



k<k* 



a 



Remark 5. Notice that if oj is rotationally symmetric, then a = 0. Indeed, since 
Xn is rotationally symmetric whenever E n is simple, this follows from (|12p . 

Remark 6. We point out that 3r(£) is a symmetric and positive operator for £ 
real (see, e.g., [TO]). We can thus generically expect that for any E > Ei fixed 
there exists e* > small enough such that HY(9) does not have discrete spectrum 
in the interval [E\, E] for any < |e| < £*; more precisely, for any S > the set 



is empty or it consists of finitely many points with imaginary part strictly negative 
(see Fig. U]). As a result it follows that the embedded eigenvalue E 6 S + of the 
untwisted model turns into a resonance when an appropriate twisting is applied. 

Concerning the assumption on the multiplicity of E we note that it is closely 
related to the multiplicity of E n , the eigenvalues of — A^,. In general, and especially 
for cross-sections with a high degree of symmetry, some eigenvalues of — A^, might 
be higly degenerated. However, this degeneracy is unstable under small pertur- 
bations of the domain, 26J. In fact, if the boundary of u> can be C k — smoothly 
embedded in M 2 with k > 4, then the eigenvalues E n of — A^, are generically simple, 
see [IS]. 

Remark 7. In jT?l Thm. 1] Grushin has obtained a similar asymptotic expansion 
for the ground state in mildly twisted tubes with transversal potential. He then 
proved, under certain assumptions on to, the absence of discrete eigenvalues for e 
small enough. We would like to mention that there is one important difference 
between our result and that of [14], [17]. Assume that the boundary of uu is 
sufficiently regular (e.g. C 1 — smooth) and replace the Dirichlet boundary condition 
at dQ. by the Neumann one, let us denote the resulting operator by HY ,N Since V 
depends only on x, the eigenfunction associated with the lowest eigenvalue of H ' N 
will be given by ipi (x, s) = ctpi (x), where c is a constant. This follows from the fact 
that Xi is constant in the case of Neumann boundary conditions. Consequently 
we have d T ipi = and therefore, using ipi as a test function, we find out that 



This means that, contrary to the Dirichlet case, the lowest eigenvalue of H ' will 
not be removed by the twisting, even for very small V . On the other hand, the 
above analysis shows that the effect on the embedded eigenvalues will typically 
occur also in the Neumann case, since the eigenfunctions of H ' associated with 
the embedded eigenvalues are not constant in s. 



In the previous section we have seen that the embedded eigenvalues under the 
influence of twisting generically turn into resonances. However, Theorem [TJ docs 
not a priori say that the imaginary part of the resonances is strictly negative. In 




5. One concrete model 
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E l E 2 E 3 




Figure 4. The discrete spectrum of LLY(8), for e ^ small 
enough, consists of two parts; the first part is given by the real 
and simple eigenvalues (full circle) below E±, the second one is 
given by simple eigenvalues with real part larger that E\ and with 
imaginary part strictly negative (empty circle). 



this section, we will show on a concrete model that for mildly twisted waveguides 
one can guarantee the negativity of the imaginary part of a chosen resonance. 

To make this problem simpler we would like to consider a concrete model, in 
which V acts as a Dirac delta potential. However, as the Dirac delta potential is 
obviously not dilation analytic, see Assumption [Bl we will approximate it by the 
sequence 

V ^ = "o TJ( \"> V > > ( 15 ) 

2 cosh \vx) 

which converges to the delta function at zero as v — > oo in the sense of distributions. 
Moreover, to be able to give some quantitave results we assume that a(x) = x. 

Proposition 1. Let a(x) — x and assume that the embedded eigenvalue E = 
E2 + [L\ of the operator LL " is simple and that E2 — E\ > i. Here V v is given by 
(I15[) . Then in the vicinity of ' E there is an eigenvalue E{e,v) of HY"{&) with the 
imaginary part given by 

<5E{e,v) = -e 2 a{v) + 0(e 3 ) , (16) 

where 



— E — - 

hm a (u) = \C 1 \ 2 ^ C x = {xi,d TX i) L ^)- (17) 

v— >oo [E2 — E\) 

In particular, if C\ 7^ then the coefficient a{y) is strictly positive for v large 
enough. 

Remark 8. We note that with this choice of a the essential spectrum of ELY (6) 
will depend on e. This fact is not usual in complex scaling methods; however this 
inconvenience does not affect the existence of the resonances since the correction 
of the essential spectrum is small for e small and thus the non-real eigenvalues of 
LLY(9) will stay far enough from its essential spectrum. In fact, it could be possible 
to avoid this fact by choosing 

a(x) =x, Vx G [-X, X] and a(x) = , V|sc| > 2X, (18) 

for some X>1, We note that with this choice of a the essential spectrum of 
HY (9) will not depend on e and, furthermore, since the wavefunction ip\ decays 
very fast as |x| grows then, in order to compute pip and (|33|). we can practically 
take a = 1 and a = 0. 
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Remark 9. Equation (|17|) implies that if Ci 7^ 0, then the imaginary part ^SE(e, v) 
will be strictly negative for suitable values of v end e. Indeed, first we take v large 
enough so that a{v) > 0. Then we keep this v fixed and take e = e(u) small enough 
so that the remainder term 0(e 3 ), which also depends on v, becomes smaller than 
e 2 a{v). Then 

%E(e,v)<0, (19) 

which means that the twisting pushes the eigenvalue E(e, v) down in the complex 
plane, making thus the lifetime of the corresponding resonance shorter. 

Remark 10. Note that the assumption E2 — E\ > j is made only for the sake of 
simplicity. It guarantees that E is an embedded eigenvalue of Hq" for any positive 
v. If E2 — Ei < j, then we would have to consider an eigenvalue E = Ek + [J-i for 
some k large enough such that E is embedded. 

Remark 11. The coefficient C\ depends on the geometry of lo. An integration 
by parts shows that C\ = — (d T xi, X2)l 2 (uj)- If is therefore easy to see that 
C\ = for rotationally symmetric lo. However, it is not a priori guaranteed that 
C\ 7^ whenever lo is not rotationally symmetric, although we are not aware of any 
counter-example. Let us only mention that C\ 7^ for certain domains lo. Indeed, 
for lo = [0, a] x [0, b) with a > b > 0, we have 



Xi{y,z) 




A 1 1 explicit calculation then shows that 

5.1. Proof of Proposition [TJ Equation (|16p follows directly from Theorem [T] 
The rest of the proof will be given in two steps. 



Spectrum of h v = —d 2 + V v . Following [3TJ §23] we set 



1 

t = - 
2 



I W 1 + - 



The eigenvalue problem h u tf>j = Hjipj admits solutions 

~\ 2 



-(2j-l) + Wl + 



, 1 < j < t + 1 



with associated eigenfunctions 



9j 



(x) = (l-£ 2 p /2 F 



ej -t,ej + t+ 1, ej + 1, -(1 - 



where 



£ = tanh(^x), ej 



V^j _ 1 



-(2j-l) + Wl + 



(20) 



(21) 



(22) 



(23) 
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F denotes the hypergeometric function and ej — t = j — 1. In particular, when 
v ^> 1 then t ~ ^ -C 1 and the spectrum of /i consists of only one eigenvalue 

2 



Ml 



T 



-1 



1 



1 



with the associated normalized eigenvector 



<fii(x) = [l — t&nh 2 (vx] 



We recall that the absolute continuous spectrum of the operators 



(24) 



(25) 



-d x + V v [x) and h c 



-St 



where 8 denotes the Dirac's delta at x = 0, coincides with the positive real axis: 

<Jac{K) = 0- ac (hoo) = [0, +Oo) 

and that, see e.g. [H Thru. 3.2.3], h v — > /ioo as f — * +oo in the norm resolvent 
sense: 

lim 11^(0-^(011 = 0, 3C>0, 

v — >oo 

where ?v(0 = [C — h u ] 1 and roo(0 — [C — ^oo] • Furthermore, making use of 
the same arguments as in [1] §3.2], it follows that for any rapidly decreasing test 
function ip 

lim (ip, MO - roo(0] ¥>W) = 0, > 0, (26) 
uniformly for ( belonging to a compact set [a, b] x i[0, c] for any a, 6, c > 0. 
Computation of the coefficient a. For v — ► oo we have 

^(M = {mi=-|+C(^ 1 )} . (27) 

We take ^ large enough so that for the set E = {£ = = /ij + of eigenvalues 
of holds 



£ w = Ex - - + Oii;- 1 ) <E!< E ia = £ 2 - i + Oiv- 1 ) < E 2 . 

Then we apply the perturbative theory to the embedded eigenvalue 

E = Ei t 2 = E% + fix 
with the associated eigenvector 

ip(x,y,z) = ipx(x)x2(y,z) . 

In such a case k* = 2 and 



(28) 
(29) 
(30) 
(31) 



lim 3 I ^2 (dk,r„(E - E k - ip)d k ) L 2 {R) - (ip, UY"ip) L 2, n (ip) 1 



fe=l,2 



where rv(0 is the resolvent of h v and 

d k(x) = (Xk,U^i(j) L 2 {Lu) 



(32) 
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An integration by parts shows that (x2, d T X2) l 2 (w) = 0. We thus get 

~e 2 ((p 1 ,a 2 (p 1 ) L 2 {R) (x2,d^.X2)L^(u J ) 

= —C E 2 , 

where 

Co = ( l Pi,a 2 <Pi)L2(R)(X2,d 2 X2)L2(u J ) = (X2,d 2 X2)mu J ) ■ 
Furthermore, from (f3"2")) and we obtain that 

di(x) = -2eCi d x f! - e 2 C 2 ipi, d 2 {x) = ~e 2 C cpi (33) 

where 

Ci = (xi,9 t X2)l 2 {u}) and C 2 = (xi, 9^X2} ■ 

Collecting all these facts and keeping in mind that E = E 2 + p\ and <pi is the eigen- 
function of h v with eigenvalue p\ we get after some tedious, but straightforward 
calculations, that 



hm + 9 iy{d 2 ,r v (E- E 2 - ip)d 2 ) L 2 {R) - \(^,U^ip)\ 2 L2{n) (ip) x | = 0. 
This implies 



3a 2 = - lim 3 (d\,r v (E ~ Ei - ip)di) L 7.(v ) (34) 
p— >o+ 

: - lirn %[4e 2 \C 1 \ 2 {d x <p 1 ,r u {E-E 1 -ip)d x i Pl }] +0(e 3 ). (35) 
p— >o+ 

We now pass to the limit v — » oo which implies 

Mi-^-^) Pi -* = e~ N/2 , hv-^hoo, (36) 

where the last limit is reached in the norm resolvent sense. We recall also that the 
resolvent [C — /loo]^ 1 has the kernel given by 

JCq{x 1 x') — 1C®(x, x') + ICAx, x') , 

where 

K°(^0 = ^l*-*'l, K\(x,x>) = -l-J_. e ^W^'l], C = fc 2 ,^>0, 

see pp. In our case 

C = k 2 = E-E 1 -ip = E 2 -E 1 +p 1 -ip. (37) 

We will denote by and K.^ the integral operators with the kernels }C®(x, x') and 
K,^(x,x') respectively. Note that d x ipi(x) is an odd function, which implies that 
K.^d x Lpi = 0. Now hi, ► /ioo as in (|26[) . Since E E\ is not an eigenvalue of \i v 
for any v large enough (in fact E — E± belongs to the absolute continuous spectrum 
of the operators h v and hoc) and c?i is an exponentially decreasing function as 
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\x\ — > oo, we can pass to the limit v — ► oo replacing r v (E — E\ — ip) on the r.h.s. of 
flS} by KP C and <pi by <j>: 

lim 9a2 = — lim lim 3 (d%, r v (E — E\ — ip) d\)r,2 tm (38) 
- - lim Of Ue 2 \C 1 \ 2 {d x 4>,roo(E-E 1 -ip)d x 4>)} +0(e 3 ), 

p— +0+ 

where the remainder term is uniform with respect to p. An explicit computation 
then gives 



^' C * W [1+4(^2-^+/!!)] 2 



In view of (|36| and (f39|) we get 



E2 ~~ E\ — j 

lim 3a 2 = -£ 2 |Ci| 2 V 



(£ 2 -£i) 2 
The proof is complete. 
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